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ABSTRACT 

Gamma ray bursts are often modelled as jet-like outflows directed towards the ob¬ 
server; the cone angle of the jet is then commonly inferred from the time at which 
there is a steepening in the power-law decay of the afterglow. We consider an alter¬ 
native model in which the jet has a beam pattern where the luminosity per unit solid 
angle (and perhaps also the initial Lorentz factor) decreases smoothly away from the 
axis, rather than having a well-defined cone angle within which the flow is uniform. We 
show that the break in the afterglow light curve then occurs at a time that depends 
on the viewing angle. Instead of implying a range of intrinsically different jets - some 
very narrow, and others with similar power spread over a wider cone - the data on 
afterglow breaks could be consistent with a standardized jet, viewed from different 
angles. We discuss the implication of this model for the luminosity function. 
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1 INTRODUCTION 

There are strong reasons for suspecting that the emitting 
plasma of y-ray bursts (GRBs) is geometrically beamed in 
a cone. The energy requirements can then be reduced be¬ 
low the exorbitant levels that isotropic emission would im¬ 
ply (Kulkarni et al. 1999) and in most models for the long 
bursts it is in any case a natural expectation - borne out 
by simulations (MacFadyen & Woosley 1999; MacFadyen, 
Woosley & Heger 2001) - that the relativistic outflow from 
a central engine should be collimated along a channel that 
opens up along the rotation axis of a massive star. Although 
the relativistic MHD that gives rise to jets is uncertain, and 
surely very complicated, most discussions of the radiation 
from gamma ray bursts (and their afterglows) has postu¬ 
lated a jet with a well-defined angle (Meszaros & Rees 1997), 
though this angle may differ for different bursts (see, how¬ 
ever, Meszaros, Rees & Wijers 1998; Salmonson 2001). We 
discuss here an alternative model where the jet, rather than 
having a uniform profile out to some definite cone angle, has 
a “beam pattern” where the power per unit solid angle (and 
perhaps also the initial Lorentz factor) is maximal along 
the axis, but drops oft gradually away from the axis. This 
would be expected if there is mixing and entrainment from 
the borders of the funnel. We discuss the expected time- 
dependence of the afterglow if it is triggered by a jet with 
this more general profile. We conclude that the time of the 
observed break (which, for a uniform jet viewed along its 
axis, depends on the cone angle; see, e.g., Rhoads 1997) in¬ 


stead depends on the angle between the line of sight and the 
symmetry axis: such a jet viewed nearly head-on simulates a 
narrow uniform jet, whereas the afterglow from the same jet 
viewed more obliquely would simulate a wider uniform jet. 
GRB have been proposed recently to be explosions releasing 
a standard power that can be either injected in very different 
jet opening angles or distributed within the jet in some uni¬ 
versal emission diagram, (Postnov et al. 2001). While Frail 
et al. (2001, hereafter FOl), support with their data the 
first interpretation, we show that their observational results 
could instead be attributed to a more standard set of objects 
viewed at different angles to their symmetry axis. 


2 DYNAMICS AND BREAK TIME 


We suppose that all long GRBs have jets with a standard 
opening angle 6j, total kinetic energy and beam profile for 
0 < 9 < Oj. We consider a relativistic outflow where both 
the bulk Lorentz factor and the energy per unit solid angle 
depend as power law^ on the angular distance from the 
center 6 



0 <9 < 9 ^ 
9c<9 <9j 


( 1 ) 


^ We concentrate here on an power index —2 motivated by FOl 
correlation. This is not, however, a necessary ingredient of the 
model as discussed in §6 (see also Fig. 3 and 4) 
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and 

r^ f r. o<e<e. 

1 ar>0 e„<e<ej, 

where 9c is introduced just for formal reasons to avoid a 
divergence at 0 = 0 , but can be taken to be smaller than 
any other angle of interest. A lower limit to this angle is 
9c > l/Fmax ~ 10 “® degrees, where Fmax ~ 10 ® is the max¬ 
imum value to which the fireball can be accelerated (Pi- 
ran 1999). The power law index of F, ar, is not impor¬ 
tant for the dynamics of the fireball and the computation 
of the light curve as long as F(t = 0, 0) = Fo(0) > 9~^ and 
Fo( 6 ') S> 1, V0. Nevertheless, it plays a role when we want 
to calculate the fraction of GRBs-afterglow without prompt 
7 -ray emission, as we discuss in § 6 . 

Consider an observer at an angle 9o < 9j with respect 
to the axis of the jet. He measures an isotropic equivalent 
energy Eiso ~ 4ne{9o) from the 7 -ray fiuence. If also the 
afterglow emission is dominated by the component pointing 
the earth, he will infer 9o as the half-opening angle of the jet 
by means of the break time in the light curve, tt, (Sari, Piran 
& Halpern 1999) 9j oc (Eiso/n)~^''^ , where n is the ex¬ 
ternal medium density. Since the total energy inferred from 
all viewing angles is Etot — 211 9^t =const, the observer will 
derive the same conclusions obtained by FOl and Panaitescu 
& Kumar (2001, hereafter PKOl) of GRBs as fireballs with 
the same total kinetic energy but very different jet apertures. 

To evaluate how the contributions of the other compo¬ 
nents add to the light curve from the zone with 9 9o we 
calculate when their beamed emission include the observer 
direction, (when 9 — 9c, < ^) and which energy per unit 
solid angle they have at that time compared to e{9o) = to- 
We show that under the assumptions of Eq. the afterglow 
light curve is indeed dominated by the fireball element along 
the line of sight: neither the “core” of the jet with 9 9o 
nor the regions with 9 9o make substantial contributions. 

For an effective though simplified discussion we con¬ 
sider only three components of the jet corresponding to 
9 = 9i 9o, 9 = 9o and 9 — 9^ ^ 9o- We call them 
“cone i”, where i — l,o, 3, respectively. This approach is 
justified by the fact that only the very inner parts of the jet 
with e ^ Co and the much wider (^ 6 ^^ tt^o) outer parts 
could contribute and substantially modify the light curve 
of cone o. We approximate cone 1 as a relativistic source 
moving at an angle 9o with respect the observer, while the 
observer is approximately considered on the symmetry axis 
for cones o and 3. The cones, that are not causally con¬ 
nected (Fo > 9~^ ), evolve independently and adiabatically 
in a constant density medium and spread sideways when Fi 
drops below 9~^ (Rhoads 1997). We consider relativistic lat¬ 
eral expansion of the cones so that their geometrical angles 
9i grow as Si ~ 0i(t = 0) -I- Since the dynamics of the 
fireball components changes at tf,, also F decreases with time 
differently before and after the break time. For t <tb 


I Fo(^)-®/" for Si 
i Fo(j^)-®/® for 9c and S 3 , 


( 3 ) 


where td is the deceleration time at which F = Fo/2 and the 
jet begins to decelerate significantly. From Eq. ^ to Eq. ^ we 
drop the subscript i for an easier reading. For t > tb 



This on-axis calculation is valid for all Si, because Fi ^ 1/So 
soon after the break time and the emitting plasma enters the 
line of sight: from this moment we can consider the observer 
to be on the cone axis. The lateral expansion starts at a time 


r 9^/^Vy^td for Si 
\ S*/®F®/®td forSoandSs, 


( 5 ) 


where the deceleration time, td, depends on 


r Fq ^^®(1 —/3cosSo) for Si 

^ Fq for So and S 3 , 


( 6 ) 


and consequently 
Rb oc 

Therefore from Eq. ^ it follows that J? 6 =constant and 
/ S^/® 9l oc 9l for Si 

oc I ^8/3 ^1/3 ^ 


( 7 ) 

( 8 ) 


Cones o and 3 are hollow but they spread only outwards 
because the inner components which have already expanded 
(for cone 3) or are about to do so (for cone o) have higher 
pressure. Eq. ^ shows also that when Fi ^ S^^ and the 
part of the blast generated by the core of the jet becomes 
visible to the observer, cone o has just started spreading. 
The contribution of region 1 to the light curve of the cone 
with So is not dominant because when Fi ^ 9 c^ the energy 
per unit solid angle is comparable to to 



The regions with 9 ^ 9o spread at later time, so when 
F 3 ~ 9^^ the energy per solid angle of cone o has already 
decreased with time but it remains always comparable to £3 



where tb- is the break time for cones i and we used Eq. ^ 
and Eq. 

Therefore the regions with 9 9 o, (or those with S 

substantially above So), don’t determine the overall shape 
of the light curve, because they are hidden by the dominant 
emission from the component along the line of sight. In prin¬ 
ciple, the superposition of many components of energy given 
by Eq. ^ and may give rise to a sizable contribution. As 
shown by the more detailed calculation of § 3.2 (see Eig 2), 
this is not the case and the only effect is to delay by a factor 
of 2 the time of the break. Gonsequently the energy per unit 
solid angle that is measured modeling the afterglow is to and 
the time break depends only on the viewing angle So. 


3 LIGHT CURVE CALCULATION 

In order to compute a reasonably accurate lightcurve in our 
model, we adopt the following approach. We divide the in¬ 
homogeneous jet in N hollow cones (all but the very central 
one), each characterized by energy per unit solid angle and 
Lorentz factor given by Eq. We compute an approximate 
lightcurve for each sub-jet using three asymptotic behaviors. 
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t (days) 


Figure 1. Light curves of an homogeneous jet (with uniform F 
and e) as observed from different viewing angles, the geometric 
opening angle of the jet is thetaj = 1°. From top to bottom, with 
different line styles, off-axis angles 9o = 0, 2, 4 and 8 are shown. 

For the cones with 9 6o and 9 9o 'fie adopt the 9 9o 
and 9 ^ 9o approximations described in § 2,. For the cones 
defined by — 1/F < 9 < 9c, + 1/T which point towards the 
observer since the beginning, we again consider the observer 
along the symmetry axis, but the emission is calculated in a 
filled cone geometry, instead of the hollow one adopted for 
the other cases. For 9 > 9o — 1/T we use the usual after¬ 
glow theory to perform the lightcurve calculations, while for 
9 9c, we generalize it for an off-axis observer as explained 

in the following section. 

3.1 Radiation from an off-axis homogeneous 
fireball 

Consider a uniform jet with Lorentz factor F and initial half¬ 
aperture 9j. The radiative process is synchrotron emission 
(Meszaros & Rees 1997) and we concentrate on the power 
law branch of the spectrum between the peak, i/m, and the 
cooling, j/c, frequencies. The observed flux at a frequency v, 
time t and viewing angle 9o is 

F{v,t,9o) oi Ac 1'5^{T,9o), (11) 

where is the emitting area and 5 = (F (1 — /3cos6))“^ is 
the relativistic Doppler factor. I' is the comoving intensity 
at the comoving frequency v' = v/S and at the comoving 
time t' = St 

I' = I'{v'cr„t')(^\ “ oc (12) 

\ / 

Due to the relativistic beaming, the observed flux depends 
on the observer angle. In particular, if 6^0 ~ 0 


Figure 2. The light curves of an inhomogeneous jet ob¬ 
served from different angles. From the top we show do = 
0.5, 1, 2, 4, 8, 16°. The break time is related only to the ob¬ 
server angle: cx. d/,. The dashed line is the on-axis light curve 

of an homogeneous jet with an opening angle 29o and an en¬ 
ergy per unit solid angle e(6o). The blow up is the time range 
between 4 hours and 1 month, where most of the optical obser¬ 
vations are performed. Comparing the solid and the dashed lines 
for a fixed 9o, it is apparent that we can hardly distinguish the 
two models by fitting the afterglow data. The slight flattening of 
the lightcurves just before the break is likely to be due to the 
numerical approximations we adopted (see § 3). 


In this case the jet is initially a source moving at an angle 
9o with the observer. At the break time F ~ 9~^ and then it 
starts decreasing very fast (Eq. until 9o — 9j < f and we 
can again use the approximation that the jet is viewed on- 
axis. So the intensity, the change in the slope and the break 
time of the afterglow depend on the observer viewing angle. 
In Fig. 1^ we show the resulting off-axis lightcurve for a jet 
with opening angle 9j = 1° as viewed with off-axis angles 
9o = 0, 2, 4 and 8 degrees. The flat part of the lightcurves 
corresponds to the time interval between the beginning of 
the jet spreading and the time in which the jet enters the 
line of sight. In this time interval our approximations are no 
longer valid and we simply model the lightcurve with a flat 
component connecting earlier and later times. This crude 
approximations is likely responsible for the slight flattening 
of the final lightcurves (see below and Fig. just before the 
break. 


3.2 Inhomogeneous jet 

Using Eq. and Eq. ^ we can compute the total light 
curve as a sum over the fireball components according to 
the scheme described at the beginning of this section. 


F = Fon — rr 
while for 9o > 9j 


7'(2^)^ vt 


F = Foff ~ 


ti{R 9) I {Y(i-pcoseo)'^ t < tb 
H^2j/(2r)3 t>tb. 


4 #) 


(13) 


(14) 


L{v,t,9o) ^e<eo-Y, 


Foffi 




+ 

Fon Fon 


(Sira)] , 


( 15 ) 


where 9in is the inner edge of the cone and Fon{9in) is the 
light curve for a jet with e{9), T{9) and half-opening angle 
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6in. N is the number of jet components that are beamed 
around the line of sight at small times. 

The results of this calculation are shown in Fig. ^ for 
different viewing angles. On top of the solid lines we plot 
with a dashed line the lightcurve of a homogeneous fireball 
with £ = Co and 9j = Bo- This shows that, for Bo < Bj, 
the lightcurve of an inhomogeneous jet can be successfully 
modelled with the lightcurve expected for a homogeneous 
jet: F oc for t < tf, and F oc for t > tf,. 

One of the main simplifications that we made in Eq. 
is to assume that the observer is on the jet axis of the cones 
with B > Bo- The main consequence is to predict a transition 
between the two power-law branches sharper than what we 
expect from the exact integration. This is due to the fact 
that, even in a uniform jet, off-axis observers sees smoother 
transitions (see Fig. 4 of Ghisellini & Lazzati 1999). 


4 TIME BREAK-ENERGY RELATION 

Recently an important observational result on the ener¬ 
getic content of GRBjs was published by PKOl and FOl. 
They found an anti-correlation between the isotropic equiv¬ 
alent energy, Eiao = dvre and the break time in the after¬ 
glow lightcurves. FOl in particular derived Eiao from gamma 
ray fluences, Ey oc Eiao and their data are consistent with 
Ej oc They explain it in the framework of collimated, 
uniform, lateral spreading jets interacting with a constant, 
low density (0.1 cm“®), external medium. The observer is 
postulated to be along the jet axis. They assume that the 
afterglow emission steepens because at that time the Lorentz 
factor has dropped to F ~ 0 “^, so that the on-axis observer 
sees the edge of the jet and the lateral jet spreading be¬ 
comes important. In this case Eiao oc is predicted, which 
matches observation. They convert the observed break times 
in jet opening angles through the formulation of Sari et al. 
(1999). The 7 -ray energy measured and corrected for the in¬ 
ferred geometry of the jet is clustered around 5 x 10®° erg. 
They concluded that GRBs central engines release the same 
amount of energy through jets with very different opening 
angles. In this framework, then, the wide distributions in 
kinetic energy per unit solid angle, which spans 3 orders of 
magnitude, is due only to the distribution of jet solid an¬ 
gles. PKOl obtained the same results, modeling a subset 
of multiwavelength afterglows from which they could assign 
an external density, a time break and an equivalent isotropic 
energy for the fireball. Both FOl and PKOl found geometric 
angles that span an order of magnitude but strongly con¬ 
centrate around 2° — 4°. In the framework of our model 
an alternative explanation of the observed relation between 
Eiao and tj can be found. As discussed in §2 and shown in 
Fig. we can infer from observations only the properties 
of the cone pointing towards the observer and not those of 
the whole jet. So Eiao = dvreo and using Eq. ^ and Eq. we 
obtain 

Eiao oc Bo^ oc . (16) 

Then, the observed distribution of Eiao and its relation with 
th is due to an inhomogeneous jet and the possibility to view 
it from different angles. 


5 LUMINOSITY FUNCTION 

Under our assumptions, each 7 -ray luminosity corresponds 
to a particular viewing angle. The probability to see a jet 
between B and B F dB is given by 

P{B)dB oc sin BdB 0 < 0 < (17) 

(0) ~ 0.7 and the highest probability is for B = Bj. Therefore 
it is highly improbable to see a jet on axis. Consequently we 
expect more faint GRBs then very luminous ones according 
to a luminosity function 

P{y) cx 10-", (18) 

where y — log e. 

Since there are only a small number of GRBs with ob¬ 
served red-shift, the comparison of our predicted luminos¬ 
ity function with data is far from being definitive. Recently 
Bloom at al. (2001) published an histogram of the bolomet- 
ric k-corrected prompt energies for 17 GRBs. The distribu¬ 
tion is roughly flat from 6 x 10 ®^ to 2 x 10 ®“^ erg but as 
the authors emphasize this analysis applies only to observed 
GRBs with redshifts and several observational biases ob¬ 
scure the true underlying energy distribution. The main bias 
that overcasts faint GRBs is the detection threshold of the 
instruments: this sample is thus flux and not volume limited. 
Moreover redshift determination encounters more problem 
for faint GRBs. Based on (U/Una!i:)-hardness correlation, 
Schmidt (2001) derived a luminosity function without using 
any redshift. They had to assume how the comoving GRBs 
rate varies with redshift and they based their calculations 
on three star forming rate models. In this case a power- 
law luminosity function was derived, but flatter than the 
n{L) oc L~^ predicted in the simplest version of our model. 
A luminosity function not based on assumed burst rate evo¬ 
lution can be derived by measuring the burst distance scale 
through the recently discovered variability-luminosity rela¬ 
tion (Fenimore & Ramirez-Ruiz 2000, Reichart et al. 2001). 
A cumulative analysis of a sample of 220 bursts (Fenimore 

6 Ramirez-Ruiz 2000) yielded a power-law luminosity func¬ 
tion n{L) oc which compare more favorably with our 

prediction. More recently, Lloyd-Ronning, Fryer & Ramirez- 
Ruiz (2001) find that the typical luminosity of GRBs evolves 
with redshift. As a consequence, a flatter power-law index 
n{L) oc L~^'^ is obtained. We stress, however, that our de¬ 
duced luminosity function can be altered by distributions in 
total energy or geometric angles or by luminosity evolution 
with redshift. Given the somewhat contradictory observa¬ 
tional results discussed above, we conclude that more accu¬ 
rate spectral and fluence measures and a larger sample of 
bursts are needed for a proper comparison. 


6 DISCUSSION AND CONCLUSIONS 

We considered inhomogeneous GRBs jets with a standard 
total energy, opening angle and local energy distribution, 
e oc B~^. We show that this jet structure can reproduce 
the observed correlation between isotropic energy and break- 
time. In this model both measurements depend only on the 
viewing angle because the 7 -ray fluence and the afterglow 
emission are dominated by the components of the jet point¬ 
ing towards the observer at small times. Since all cones have 
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the same total energy Ejet = 27r0^e = cost, we recover the 
results of FOl and PKOl but the constrains on the geometri¬ 
cal beaming can be relaxed and an appealing more standard 
structure for all GRBs can be adopted. The jet total energy 
can be calculated from Eq. ^ 

Erotai = 2'Kt,el (^1 -p 2 In = Ejet (l + 2 In (19) 

and compared to Ejet = 2 it 6 of-o = 2 Trec 0 f the total energy 
inferred from observation, Ejet < Exotai- To give an ex¬ 
ample, for a fireball with 6 c — 1 ° and 6 j = 20 °, we have 
Erotai/Ejet — 6 , i.e. the true energy of the fireball can be 
one order of magnitude larger than what inferred with the 
models of FOl and PKOl. 

In addition (§ 5) we can derive the GRBs luminosity 
function from the probability distribution of the viewing an¬ 
gle and compare it to data. The comparison is, at this stage, 
still uncertain because more accurate spectral and fluence 
measures are required to build a volume limited sample and 
confirm (or rule out) this model. 

In any fireball model considering a jet-like GRB struc¬ 
ture, a certain number of afterglows without 7 -ray emission 
(orphan afterglows) is expected. For an homogeneous jet, or¬ 
phan afterglows are possible only for viewing angles greater 
than dj -I- 1/F. In our jet configuration a fraction of the total 
area could have a Lorentz factor lower than the minimum 
F necessary for 7 -ray radiation. Gonsequently, considering 
the same opening angle, an inhomogeneous jet could pro¬ 
duce an higher fraction of orphan afterglows than an ho¬ 
mogeneous one. This intrinsic fraction depends above all on 
the F distribution within the jet (ar in Eq and on the 
minimum Lorentz factors to produce 7 -ray radiation and 
afterglow emission. The observed number of orphan after¬ 
glows depends also on flux detection limits, GRB explosion 
rates with redshift and cosmology. An accurate calculation 
of the expected orphan afterglow rates is therefore beyond 
the scope of this paper. 

We emphasize that we implicitly assumed the radiation 
efficiency of the fireball to be weakly dependent on F or e. 
This is a plausible assumption in internal-shocks scenario. If 
it was not the case and the efficiency grew with e, a different 
relation between e and 9 should be postulated in order to 
reproduced observations. In this paper we concentrated for 
simplicity on a beam profile e oc d~^ which is consistent 
observational results but it is interesting to briefly discuss 
other power-law relations e oc 0““' (see Fig. ^). A decay 
flatter then 2 would cause two breaks in the light curve: the 
first due to the cone pointing the observer when r( 6 'o) ~ 
9c ^ and the second, at later times, when the observer sees 
the edge of the jet and r(0j) ~ 0“^. The power law index 
after the first break is flatter then because the cones 
with 60 < 9 < 9j enter the line of sight with e > eo and 
substantially modify the light curve shape. With a steeper 
decay in the distribution of e the time break and the emission 
after that break would be dominated by the jet along the 
axis rather then by the very much weaker part directed to 
the observer. The jet break would then be preceded by a 
prominent flattening in the lightcurve, especially for Oe > 
3, difficult to reconcile with observations. In this case we 
would have 7 -ray emission only for very small angles and 
the number of orphan afterglows would be much greater then 
expected from the Qe = 2 model. From Eq. ^ we can derive 



Figure 3. The lightcurves of inhomogeneous jets with different 
indices (e oc 0““'). The lightcurves have the same Co but they 
are plotted shifted by factors of 10 for clarity (they would be 
indistinguishable at small time). The arrows highlight the location 
of the two breaks for 0 < o < 2. 
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Figure 4. The relation between the indices at {Etgc oc “*) and 
Oe (e oc 0““®) from an inhomogeneous jet. The shaded regions 
show the best fit to FOl data at = 0.9 ± 0.22. 


the relation between the index and at (where Etso oc t“*). 
We obtain at = a^ 12 for Qe > 2 and at = 3 aE /(8 — aj for 
Oe < 2, when the first break is considered. In Fig. M we 
plot this relation overlaid on the interval in at allowed by 
observations (we used FOl data). We derive 1.5 ^ cte 2.2, 
at the Icr level. Further 7 -ray and afterglow observations will 
allow to constrain this parameter much better in the future. 

Besides the luminosity function discussed in § 5, there 
are several ways in which this model can be proved or dis¬ 
proved. First, we have shown that the real total energy of 
the fireball can easily be an order of magnitude larger than 
what estimated by PKOl and FOl. In this case, after the fire¬ 
ball has slowed down to mild-relativistic and sub-relativistic 
speed, radio calorimetry (Frail, Waxman & Kulkarni 2000) 
should allow us to detect the excess energy. In addition, in 
this model we naturally predict that the more luminous part 
of the fireball have higher Lorentz factors. This may help 
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explaining the detected luminosity-variability (Fenimore & 
Ramirez-Ruiz 2000) and luminosity-lag (Norris, Marani & 
Bonnell 2000) correlations ( Salmonson 2000, Kobayashy, 
Ryde & MacFadyen 2001; Ramirez-Ruiz & Lloyd-Ronning 
2002). Another constrain is given by polarization. Since fire¬ 
ball anisotropy is a basic ingredient of this model, inducing 
polarization (Ghisellini & Lazzati 1999; see also Sari 1999 
and Gruzinov & Waxman 1999). The time evolution of the 
polarized fraction and of the position angle are however dif¬ 
ferent from a uniform jet (Rossi et al. in preparation). Fi¬ 
nally, the properties of the bursts should not depend on the 
location of the progenitor in the host galaxy, and therefore 
this model can accommodate the marginal detection of an 
i?iso-offset relation (Ramirez-Ruiz, Lazzati & Blain 2001) 
only if a distribution of Oj is considered. 


ACKNOWLEDGMENTS 

We thank G. Ghisellini, J. Granot, P. Kumar, A. Panaitescu 
and E. Ramirez-Ruiz for many stimulating discussions. ER 
and DL thank the IAS, Princeton, for the kind hospitality 
during the final part of the preparation of this work. ER 
thanks the Isaac Newton and PPARG for financial support. 


REFERENCES 


Bloom J. S., Frail D. A., Sari R., 2001, AJ, 121, 2879 
Fenimore E. E. . Ramirez-Ruiz E., 2000, ApJ submitted (astro- 
ph/00041^ 

Frail D. A., Waxman E., Kulkarni S. R., 2000, ApJ, 537, 191 
Frail D. A., et al., 2001, ApJ, 562, L55 
Ghisellini G., Lazzati D., 1999, MNRAS, 309, L7 
Gruzinov A., Waxman E., 1999, ApJ, 511, 852 
Kohavashi S.. R vHe F., MacFadyen A., ApJ submitted (astro- 
ph/011008^ 

Kulkarni S. R., et ah, 1999, Nature, 398, 389 

TJovd-Fonni ng N. M.. Frver C . L., Ramirez-Ruiz E., 2001, ApJ, 


in press (astro-ph/0108200) 

MacFadyen A. L, Woosley S. E., 1999, ApJ, 524, 262 

MacFadyen A. L, Woosley S. E., Heger A., 2001, ApJ, 550, 410 

Meszaros P., Rees M. J., 1997, ApJ, 476, 232 

Meszaros P., Rees M. J., Wijers R. A. M. J., 1998, ApJ, 499, 301 

Norris J. P., Marani G. F., Bonnell J. T., 2000, ApJ, 534, 248 

Panaitescu A., Kumar P., 2001, 560, L49 

Piran T., 1999, Physics Rep., 314, 575 

Postnov K. A., Prokhorov M. E., Lipunov V. M., 2001, Astronomy 
Report, 45, 236 

Ramirez-Ruiz E., Lazzati D., Blain A. W., 2002, ApJ, 565, L9 
Reichart D. E., Lamb D. Q., Fenimore E. E., Ramirez-Ruiz E., 
Cline T. L., Hurley K., 2001, ApJ, 552, 57 
Rhoads J. E., 1997, ApJ, 487, LI 
Salmonson J. D., 2000, ApJ, 544, L115 
Salmonson J. D., 2001, ApJ, 546, L29 
Sari R., 1999, ApJ, 524, L43 

Sari R., Piran T., Halpern J. P., 1999, ApJ, 519, L17 
Schmidt M., 2001, ApJ, 552, 36 











